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Abstract. We present 1 an exact normalisable zero-energy chiral fermion 
solution for abelian BPS dipoles. For a single dipole, this solution is con- 
tained within the high temperature limit of the SU(2) caloron with non- 
trivial holonomy. 

^ , 1. The Dirac Monopole 

00 ' 

' A convenient representation for the Dirac monopole [1] is given by 

(N \ 

A = gn A V log w(x), w(x) = \x\ + x ■ n, (1) 

where w(x) is positive, but vanishes along the Dirac string pointing along 
' ■^ ■ — h, as seen from the monopole. The magnetic field is 

B = V A A = -gV (n ■ V log w(x)) + ghA log w(x) = B rcg + S str , (2) 

and using that h ■ Vlogw(x) = l/\x\ (the derivative along the direction 
of the Dirac string) is independent of n, the first term gives rise to the 
radial magnetic field associated with a magnetic point charge, whereas the 
second term represents the Dirac string, which takes care of the return flux. 
This follows from the fact that logii;(x) is harmonic, except where w(x) 
vanishes. To be specific, choosing for convenience h = e% = (0, 0, 1), we find 
B = gx/\x\ 3 + 4Tige 3 5{x)5(y)6(-z) and V • B vcg = -V • B stl = 4vr^ 3 (x). It 
gives the appropriate magnetic point charge for B TCg , but when including 
the return flux V • B = 0, as it should. 

The function logw(x) can be viewed as a potential, although from the 
point of view of the Maxwell equations it is more natural to consider 
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A) = — gn • Vlogio(x) as such. Not only S reg = VAj, but also as the 
time component of the Euclidean vector potential this choice of Aq gives 
rise to a self-dual configuration, with E = VAq = B TCg . The usefulness 
of w(x) becomes clear when one considers the massless (Euclidean) Dirac 
equation in such a background. As usual we split this into positive and 
negative chirality Weyl equations, 

H+ = D = D^ = -o»D^ H_ = D = u^D^ (3) 

where = <9 M + ieA^ is the covariant derivative and aj = iTj, whereas <7o 
is the 2x2 identity matrix. For n = e 3 a solution of D^> = is given by 

*(*) = ("^ +^ 2 ) \ogw{x). (4) 

If so desired a (spin-)rotation allows one to obtain the solution for arbitrary 
n, but to keep things simple we stick to h = e 3 , such that 

( d 3 + id di- id 2 \ ( -d 3 -di + id 2 \ 

(5) 

Using the Dirac quantisation condition eg = ±, one easily verifies 

-tfwo = V^(_ A ^)- (6) 

Since y/w vanishes along the Dirac string, we find that y/w/S.\ogw = 
(as a distribution), hence D^f = 0. Likewise |$| 2 = wVlogw ■ Vlogw = 
Aw—wA log w. Thus |^| 2 = Aw = 2/\x\ has an integrable singularity at the 
origin, and the Dirac string is invisible, as it should. Nevertheless, ^ does 
not decay sufficiently fast to be normalisable. Note that the zero-mode is 
time independent. Putting do to zero, H is precisely the Dirac Hamiltonian, 
with Aq playing the role of a Higgs field. A non-zero asymptotic value of 
Aq would lead to a mass scale and exponentially decaying wave functions 2 . 

It should not come as a surprise that existence of zero-energy solutions 
is sensitive to the sign of the electron charge (relative to g). With eg = — \ 
we have 



-iH, 



1 / d 3 \ f J\ _ 1 /Alogw(x) 



^)\d 1+ id 2 ) l °^ )■ ^ 



However, here the singularity of the Dirac string is no longer nullified but 
enhanced. Nevertheless, it can be turned into a proper zero-energy solution, 

2 A constant Ao in the Euclidean Weyl equation can also be identified with a non-zero 
chemical potential. 
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identical to ^>'a2, by replacing w{x) with l/w(x), but this has the same 
effect as changing g to —g, explaining why the new zero-energy solution is 
the charge conjugate of ^ in Eq. (4). Negative chirality zero-energy solu- 
tions cannot appear because the self-duality of implies that DD = —D 2 . 
Therefore such a solution would satisfy = 0, which is ruled out. 

It is well known that the 't Hooft-Polyakov monopole [2] allows for a 
normalisable chiral zero-energy solution of the Dirac equation [3] . The size 
of the core of these non-abelian monopoles is determined by the mass scale 
set by the asymptotic value of the Higgs field. When the core size shrinks 
to zero, so does the support for the zero-mode. Adding to the Higgs field 
(J> = ±<E> a cr a a constant element in U(l), <1? 2 = <5 — 2niz, as it appears in 
Nahm's work [4], the zero-mode remains normalisable for a finite range of z 
determined by the Callias index theorem [5]. Our solution corresponds to z 
at the boundary of this range, where the zero-mode fails to be normalisable. 
This boundary value of z is defined by det(<l> 2 ) = 0, and the abelian field 
Aq given above corresponds to the isospin component of $> 2 responsible for 
this vanishing eigenvalue. It does imply the support of the zero-mode is no 
longer confined to the non-abelian core. 

There has been another context in which solutions to the Dirac equa- 
tion in the background of a monopole have appeared in the past, namely 
that of monopole- induced proton decay (the Callan-Rubakov effect [6]). 
Boundary conditions [7] for the fermions are imposed at the core of the 
monopole to describe the scattering states in the limit where the size of 
the monopole core can be neglected, so as to properly reflect the breaking 
of B-L, compatible with the chiral anomaly. There it is assumed, as for 
the Jackiw-Rebbi zero-energy solution, that the Higgs field approaches a 
non-zero constant at infinity, which through the Yukawa coupling gives a 
mass to the fermions 3 . In our case the asymptotic value of the Higgs field 
{Aq) vanishes and as we will see, the limit of zero monopole core size can 
be taken without any approximation, but at the expense of the zero-energy 
state being non-normalisable. 

2. The Abelian BPS Dipole 

To find a normalisable zero-energy solution, we have to do something about 
the asymptotic behaviour. A natural way to achieve this is to consider an 
abelian BPS (self-dual) dipole, or bipole for short. We now profit from 
having expressed the zero-mode in terms of the function w. The bipole 

3 For a non-vanishing Higgs mass, this together with the Dirac monopole field, is all 
that is left when neglecting the core of the monopole. In the Bogomol'ny limit considered 
here, identifying the Higgs field with Aq (self-duality implied by the BPS equations [8]), 
the long range component of the Higgs field modifies the Dirac equation. 
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field is generated by 

log w(x) = log (I a; I + x ■ h) — log(|x + sn\ + x ■ h + s) (8) 

where the first term represents a monopole at x = and the second term, 
with the opposite sign, an anti-monopole at x = —sn. It is convenient to 
express w(x) as 

r 1 +r 2 -s 

w(x) = , (9) 

where f\ = x and r 2 = x + sh. This shows that the Dirac strings of the 
monopole and anti-monopole partly cancel. All that we need to check is if 
with this choice of w is now normalisable. A simple computation shows 

that l^fx)! 2 = Aw(x) = — rrv, as shown in Fig. 1 (left). This is 

indeed integrable at the location of the two point charges and at infinity, 
J Sx\^{x)^ = J d 3 xAw(x) = 4tts, using that w(x) = 1 - s/\x\ +C(l/|x| 3 ). 



Figure 1. The chiral zero-mode for the bipole of charge k = 1 (left) and k = 2 (right). 
Plotted is ^(x)! 2 , the integral normalised to 1, as a function of y^x 2 +y 2 and 2 on a 
linear scale for s = .5, cutoff at respectively |*t| 2 = 5 (left) and .5 (right). 

Having found a normalisable zero-mode for one bipole, a generalisation 
to a collection of bipoles is obvious, by taking the product of w(x) for each 
such bipole. This does not affect the property that w(x) is positive, vanishes 
along the Dirac strings, and its logarithm is harmonic elsewhere. However, 
the form of the zero-mode requires all factors w to be formulated in terms 
of the same h, which means all bipoles have to point in the same direction, 
i.e. the magnetic moments of all bipoles have to be uni-directional. It is not 
clear if this is just a limitation of our simple ansatz. For multi-bipoles, all 
separated much further than each of the individual bipole sizes (s), |\l/| 2 will 
near each bipole be of the same form as for a single bipole. However, when 
two or more bipoles coincide, or equivalently when g is bigger than the 
minimal Dirac value of |^| 2 will be suppressed along the line segment 
connecting the two charges. To demonstrate this, we note that for g = j-, 
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Wk(x) = w k (x), with w(x) as given in Eq. (9). Thus, using w(x)A log w(x) = 
0, we find 

|*(aOf = A W (x? = kV-\x)A W (x) = (10) 
which integrates to Airks. The case for k = 2 is shown in Fig. 1 (right). 
3. The Caloron 

The context in which the bipole appears in a natural way is the caloron with 
non-trivial holonomy [9], in the infinite temperature limit. The periodic 
boundary conditions in the Euclidean time direction, relevant for these 
finite temperature instantons, allow for a non-trivial holonomy determined 
by the Polyakov loop, which approaches a constant value at spatial infinity, 

Poo = lim P(x), P(x) = Pexp(/ A (t,x)dt). (11) 

|:?|— >oo JO 

With Aq playing the role of a Higgs field, a non-trivial value implies that an 
SU(n) charge one caloron splits into n constituent BPS monopoles, whose 
masses are determined by the eigenvalues of the Polyakov loop 

n 

'PL = exp(27rzdiag(/xi,/i 2 ,--- ^Vi = 0- (12) 

i=i 

arranged to satisfy fi± < ■ ■ ■ fj, n < fJ, n +i = /ii + 1. The constituent masses 
8ir 2 v rn , with v m = (// m +i — fi m )/(3, add up to 87r 2 //3 such that the action 
equals that of a charge one instanton. The presence of these constituents is 
easily established from the formula [9, 10] 

TrF a 2 7 (x) = did 2 log V(z), Y>(x) = |tr(A t • • • Ai) - cos(2^), (13) 

A = _L ( rm \Pm+i\\ fcosh(2TTv m r m ) sinh(27r^ m r m )\ 
m ~ r m V r m+\ J Vsinh(27Tf m r m ) cosh(2iru m r m )J ' 

where we introduced r m = \x - y m \ and p m = y m - y m -\ (r n+ i = n and 
y n +i = yi), with y m the location of the m th constituent monopole with a 
mass 87r 2 ^ m . 

The basic ingredient in the construction of caloron solutions is the 
Greens function f x defined on the circle 4 , z 6 [0, /3 _1 ], satisfying [9, 10] 

(J~J- Z ~ *) + ^ Z ) + ^ E S ( Z ~ Vm/P)\Pm\^ fx{z, Z ) = S(Z-Z'), 

(14) 

4 For example TrF„ y (x) — —d^d^ logdet f x , leading to the result of Eq. (13). 
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where r 2 (x;z) = r 2 n for z G [p m /(3, Pm+i/P]- The variable z can be in- 
troduced through Fourier transformation with respect to time, where the 
Fourier coefficients are related to the ADHM data [11] of instantons, peri- 
odic up to a gauge rotation with Voo (giving the solution in the so-called 
algebraic gauge) . This is in one-to-one relation with the Nahm transforma- 
tion [12]. For fin/p <z'<z< fi m+ i/P Ux{z',z) = f*(z,z') for z < z') 
the explicit result [9, 13] for the Greens function can be expressed as 

7re 27rit(z-z') 

f x (z, Z) = (v m (zf)\Am-l ■ ■ ■ AlAn ■ • • An ~ e- 2mt \a 2 V m (z)) , 

r m ip 

(15) 

where the spinor v m (z) is defined by 

„ smh[2^(z - /i m /(3)r m } \ , . 

Vm{Z > ~ { cosh[27r(z - ti m /f3)r m ] J' (W) 

The chiral Dirac, or Weyl equation can be solved with the boundary con- 
dition ^> z (t + j3,x) = ex.-p(2TTizP)Voo^> z (t, x) (in addition to the two com- 
ponent spinor index, there is now also a colour index). With z = one 
obtains the finite temperature "anti-periodic" fermion zero-mode, and for 
z = the "periodic" zero-mode. 

To be specific, for the SU(2) caloron we have [12 = ~Pi = flu and 
\pi\ = \P2\ = 7fp 2 /P (p is the instanton scale parameter). The gauge field and 
zero-mode can be expressed in terms of the functions <f)~ l = l—p 2 f x (u>, to) / fl 
and x = P 2 fx(u, — u)/fl. In the algebraic gauge, choosing Voo = V%^ and 
the constituents along the z-axis (by proper combinations of gauge and 
space rotations this can always be achieved) [9] 

A a = ^a T T3<9 7 lo g + ^ Re {(vlry ~ "&y)( r l + iT 2)dfX) , (17) 

with fja~fO- a = <7[ a <J 7 ] the anti-selfdual 't Hooft tensor, and [14] 

<*> -£bJ«* ((* - 1 ) ( t : t ) /•(-■ «>) ■ 

(18) 

Particularly simple and valid for arbitrary SU(n), is the expression for the 
density of the fermion zero-modes [13, 14] 

\^ z {x)\ 2 = -^ 2 fl)^d 2 J x (z,z). (19) 

The limit fl — > can be seen as a dimensional reduction and only the 
time independent field components are expected to survive in this high 
temperature limit. It is therefore more appropriate to consider the periodic 
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gauge. For general z this periodic gauge is obtained by applying the gauge 
transformation g(t) = e~*Af ; where A™ = 2iri(u)Tz — z)S a o, such that 
* ' z = g(t)^ z an d the new gauge field A are now periodic, with 

(20) 

We find that [9] lim^o x( x ) = an d lim/j-^o <fi 1 { x ) = w(x) with w(x) as in 
Eq. (9) . The resulting abelian gauge field splits into an isospin up and down 
component, decoupled in the Weyl equation. For z = this gives a mass of 
ttui to both isospin components, which contribute equally to the density, 
and the zero-mode is supported entirely at y±, whereas for z = the 
mass is TTiy 2 , but the zero- mode is now supported entirely at y 2 . For other 
values of z the mass will depend on the isospin component, but as long as 
z / ±uj the zero-mode remains localised to either of the two constituent 
locations, jumping from one to the other when z crosses ±uj, where the 
zero-mode becomes delocalised, having support at both constituents simul- 
taneously. Indeed, for z = uj 

whereas for z = —uj, the same result follows after charge and isospin con- 
jugation, ^L w (a;) = — ^>lj(x)^a2 and ^ 2 _^{x) = ^(x)^ (with the isospin 
index made explicit). Since x vanishes in the high temperature limit, the 
only surviving isospin component of the Weyl zero-mode is the one for 
which the asymptotic value of vanishes, and for which the zero- mode 
is time independent. This non-trivial isospin component agrees (up to an 
irrelevant factor —2i\ip) with Eq. (4), for w(x) = </>~ 1 (x). 

For the high temperature limit to be smooth and unambiguous, it was 
essential that 4> be time independent, cf)~ l Alog(p = 0, and that x = 0- 
It is interesting to note that, imposing self-duality on Eq. (17), leads to a 
natural generalisation at finite temperature 

4>- 2 d 2 a log<t>+ K^i - id 2 ) X \ 2 + 1(80- id 3 ) X \ 2 = 0, (22) 
+ id 2 )<f> 2 (d 1 - id 2 ) X + <f>-\do + id 3 )<p 2 (d - id 3 ) X = o. 

In principle, but not in practise, this can be used to define (j) and x- In- 
terestingly these equations also appear when formulating self-duality in 
the so-called R-gauge introduced by Yang [15], after a suitable Backlund 
transformation [16]. 

We will end with a few words on the case of SU(n > 2). In the high 
temperature limit the zero-mode is again exponentially localised at one of 
the constituents [13]. This we can read off from Eq. (19), using the explicit 
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expression for f x (z,z) given in Eq. (15). When z passes through ii m /f3 
the zero-mode jumps from one constituent location to the other, and only 
for these values of z the zero-mode will delocalise, with the proviso that 
it will only "see" two out of the n constituents. In the periodic gauge, 
diagonalising Aq at infinity by a constant gauge rotation, we have Aq — > 
27udiag(^i//3 — z, — z, • • • , fJ, n /P ~ z ). The only non- vanishing colour 
component of the fermion zero- mode is the one for which Hj/(3 — z = 0. 
The resulting configuration is again that of the bipole in section 2. 

4. Discussion 

An interesting question is if there is some physical significance to the zero- 
modes, like for chiral symmetry breaking in the effective description of QCD 
in terms of monopoles obtained through abelian projection [17]. But first 
we would like to better understand how to go beyond the case where the 
magnetic moments of the bipoles are no longer parallel. A natural setting 
in which to address this particular question would be through the study 
of charge k > 1 calorons. It is not clear if in the high temperature limit 
a global abelian embedding for the case of non-parallel magnetic moments 
exists. Even when the magnetic moments are parallel, one would expect k 
independent zero-energy solutions, of which we have only provided one. 

What led us to the results presented in this paper, was an attempt to 
solve the Weyl equation in the background of the abelian gauge field [18] 
that is obtained from the Nahm transformation of an SU(2) charge one 
instanton on T 3 x R. This self-dual abelian gauge field is described by two 
bipoles on T 3 , but unfortunately with anti-parallel orientations. Further- 
more, two zero-modes are required in order for the (inverse) Nahm trans- 
formation to reconstruct the SU(2) instanton, for which z is to be identified 
with the time. When the zero-modes are localised at the monopole singu- 
larities, it can be shown that the resulting SU(2) gauge field is abelian. This 
describes the asymptotically flat connections (F = 0) of the instanton on 
T 3 x R, required for the integral over the action density to be finite, and 
equal to 8ir 2 . Thus, to obtain genuine non-abelian behaviour the zero-mode 
has to become delocalised for certain values of z. Here the gauge field is no 
longer flat (F / 0) and the density will have a maximum, in accordance 
with a general relationship between holonomies (with respect to each of 
the generating circles of the manifold) of a self-dual gauge field on the one 
hand and the constituent locations of the Nahm dual gauge field on the 
other hand. This dual relationship has been verified in a careful numerical 
study [19]. Similar results have been seen [20] for instantons on T 2 x R 2 . 

It is the applications to the Nahm transformation on a torus that has 
been our prime motivation for studying this problem. Taken out of this 
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context and the context of the caloron, our exact result for the chiral zero- 
mode is so simple we believe it is worthwhile to share it with the reader. 
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